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Abstract 

>  Fourier  representations  of  the  equations  of  motion  for  a  rotating 
stratified  fluid  are  obtained  both  in  the  two-dimensional  and  three- 
dimensional  wave  number  space.  It  is  shown  that  these  spectral  equations 
can  be  further  transformed  to  reveal  explicitly  the  interaction  between 
waves  and  waves  and  between  waves  and  eddies  by  decomposing  the  Fourier 
coefficients  in  terms  of  three  basis  functions,  that  correspond  to  the 
"upward”  propagating  and  the  "downward  propagating  wave  modes  and  the 
zero  frequency  horizontal  eddy  mode  in  the  linear  analysis.  In  the  case 
of  two-dimensional  rotating  stratified  fluid  motion,  the  transformation 
is  made  and  full  spectral  equations  up  to  second  moments  for  these  basis 
functions  are  derived. 
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1.  Remarks 

Rotating  stratified  wave-turbulence  interaction  is  a  highly  complex 
process,  and  one  must  continually  explore  and  examine  for  the  best 
approach  to  a  satisfactory  description  of  the  process.  In  this  report,  we 
present  a  spectral  formulation  of  the  wave-turbulence  Interaction.  Such  a 
formulation  has  the  advantage  over  the  perturbative  approach  currently 
used  in  the  literature  (e.g.  McComas  and  Bretherton,  1977)  in  that 
interaction  equations  can  be  derived  for  finite  amplitude  up  to  the 
Boussinesq  approximation,  whereas  the  perturbative  approach  by  the  very 
nature  of  the  technique  must  assume  small  amplitude  for  the  derivation. 

Our  spectral  formulation  thus  can  potentially  improve  the  description  of 
the  interaction,  especially  in  the  strong  interaction  regime. 

The  spectral  equations  of  motion  in  this  report  are  obtained  in  terms 
of  three  functions  that  represent  separately  the  ’upward*  propagating,  the 
'downward'  propagating  inertial  gravity  waves,  and  the  zero  frequency 
horizontal  eddy  current;  Equations  expressed  in  this  manner  have  the 
advantage  that  the  Interaction  between  waves  and  between  waves  and  eddies 
can  be  separated  out  and  examined  explicitly.  In  view  that  the  equations 
may  be  of  general  interest,  sufficiently  detailed  derivations  are  given. 
The  derivations  also  serve  as  a  reference  for  the  spectral  equations 
used  in  the  recent  numerical  modelling  study  of  nonlinear  internal  waves 
by  Shen  and  Holloway  (1982). 

The  interaction  equations  are  derived  below  for  a  stratified 
rotating  fluid  system  independent  of  one  horizontal  axis.  This  simpler 
two  dimensional  case  is  considered  because  it  is  most  frequently 
numerically  simulated,  and  the  full  detail  of  the  interaction 
can  be  derived.  The  more  difficult  three  dimensional  case  is  discussed 
in  a  memo  by  Holloway  which  is  attached  as  an  appendix  to  this  report. 

The  fluid  system  considered  is  unbounded,  incompressible,  uniformly 
rotating  and  uniformly  stratified  in  the  mean.  The  system  satisfies  the 
Boussinesq  approximation;  that  is,  fluid  accelerations  are  small  compared 

with  gravity  and  vertical  excursions  are  small  compared  with  density 
scale  height. 
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2.  Equations  of  motion  on  a  vertical  plane 

We  shall  briefly  state  the  equations  of  motion  from  which  the 
spectral  equations  for  the  two  dimensional  case  are  to  be  derived.  The 
relevant  quadratic  integral  invariants  for  this  case  will  also  be  noted. 

The  coordinate  system  adopted  here  is  the  usual  right-handed 
rectangular  Cartesian  coordinate  system  in  which  x  and  y  denote  the 
horizontal  axes  and  z  denotes  the  vertical  axis;  the  direction  of  z  is 
opposite  to  the  local  gravity  £  .  The  local  component  of  earth's 
rotation  f_  is  along  z.  In  this  system  the  two  dimensional  motions  are 
assumed  to  be  functions  of  x  and  z  coordinates  only.  The  full  equations 
of  motion  which  satisfy  the  Boussinesq  approximation  in  two  dimension 
are  then 


3  u 

+  u3  u  +  w5  u  -  fv  * 

-  p_1 3  P  +  vV2u 

(2.1a) 
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X  z 
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W2v 
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+  u3  p'  +  w3  p'  * 

X  z 

Rw  +  DV2p' 

(2.1e) 

where  u,v,w  are  the  velocity  components  along  x,y,z,  respectively.  P 
is  the  pressure.  pQ  is  the  mean  density,  p'  is  the  perturbation 
density  defined  by  p"  »  p  -  PQ (1  -  Rz) ,  where  R  ■  !^zPl  a  constant 
mean  density  gradient,  v  is  the  viscosity  and  D  is  the  diffusivity 
e.g.,  via  thermal  conduction. 

For  the  derivation,  it  will  be  much  simpler  to  work,  with  a  reduced 
set  of  the  above  equations.  This  is  obtained  by  eliminating  the  pressure 
through  cross  differentiating  (2.1a)  and  (2.1c)  and  by  replacing  u  and  w 


with  a  streamf unction 

t|/  via.  u  •  3  $  and  w  «  -3  . 

Z  X 

The  result  is 
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where  £  “  V2^  is  the  relative  vorticity  and  J(A,B)  ■  3  A3  B  -  3  B<J  A 

x  y  x  z 

is  the  Jacobian.  Thus  the  field  is  adequately  described  by  three 
physical  variables  <J>,  v,  and  p' . 

In  later  spectral  formulation,  we  shall  consider  only  fluid  region 
which  has  periodic  boundary  conditions.  In  such  a  region  the  net  flux 
through  the  boundary  is  zero,  and  two  quadratic  integrals  are  invariant 
provided  that  dissipation  and  diffusion  are  absent.  The  first  is  the 
conservation  of  energy,  Z^ffE  dxdz  ■  0,  where 

E  *=  (|Vt|/|2  +  v2  +  N2p'2/R2)/2  (2.3) 

and  N2  *  gR/pQ  is  the  Brunt-Vaisala  frequency  squared.  The  energy 
conservation  is  obtained  by  multiplying  the  three  equations  by  \p,  v,  and 
P'»  respectively,  summing  and  then  integrating.  The  other  invariant 
quantity  is  the  potential  vorticity 

PV  -  (O  v  +  f )  3  (p  +  p')  -  3  v3  o']  (2.4) 

X  Z  Z  X 

which  is  obtained  by  takiug  3^  and  3 z  of  (2.2b)  alternately  and  by  applying 
the  same  operation  to  (2.2c)  and  then  forming  the  product  shown  in  (2.4). 
This  is  a  special  case  of  the  Ertel’s  theorm.  These  invariants  are 
useful  as  diagnostics  for  numerical  simulation  study  and  also  suggest 
useful  second  moment  quantities  for  statistical  study  (see  eq.  3.15  and 
Sec.  3-D).  Lastly,  we  note  the  existence  of  a  third  invariant  for  the 
special  case  of  f  -  0.  In  this  case,  the  quantity  p'C  is  also  conserved. 
This  quadratic  invariant  can  be  obtained  by  multiplying  (2.2a)  and  (2.2c) 
by  p'  and  5,  respectively,  and  then  integrating. 


I 


» 


-  4  - 


3.  Derivation 

The  derivation  will  proceed  in  three  steps:  a)  Fourier  transform 
the  equations  of  motion;  b)  obtain  the  basis  functions  for  the  Fourier 
coefficients;  and  c)  construct  the  second  moment  spectral  evolution 
equations.  In  step  b,  the  basis  functions  refer  to  the  functional 
representation  of  the  two  inertial  gravity  wave  modes  -  i.e.,  the  'upward' 
and  the  'downward'  propagating  modes  -  and  the  zero  frequency  current  mode 
in  the  linear  wave  analysis.  Ue  choose  such  representation  so  that 
interactions  between  waves  and  eddies  can  be  examined  separately. 


A.  Fourier  representation 

The  domain  over  which  we  obtain  the  Fourier  expansion  is  assumed 
to  have  dimension  2v  x  2tt  .  The  choice  is  for  convenience  since  the 
dimension  can  always  be  rescaled  later  to  the  size  of  interest.  Thus 
in  terms  of  discrete  set  of  wave  vectors,  Fourier  expansions  for  iji.v.p' 
are 

(3.1a) 
(3.1b) 
(3.1c) 


<P  m  £  \  exp(ik*x) 
v  -  £  vk  exp(ik*x) 
p'm  £  p£  exp(ik*x) 


with  4^  *  ,  etc.,  since  the  physical  variables  are  real. 

Substitution  of  the  above  into  (2.2),  the  equations  for  Fourier 
coefficients  are 


W  +  vk2<J>k  +  i(fkzvk  +  gkxP'  )/k2  - 

jE  .  (k  p  -  k  p  )(q2/k2)  \(< 
£+£=U  xry  y*x  M  r£  £ 


3tvk  +  vk  v>*  +  1  fk-^  " 


z  h 

J  .  (k  p  -  k  p  )  ♦  v 
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(3.2a) 
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B.  The  Basis  Functions 

We  will  obtain  linear  solutions  to  the  equations  (3.2  ),  generalize 
the  solutions  to  basis  functions  and  then  find  the  transformation  between 
the  Fourier  coefficients  of  physical  variables  and  the  basis  functions. 

Removing  the  nonlinear  terms  from  (3.2  )  and  substituting  the 

AAA 

solution  v^,  p^}  exp(iwt)  into  the  linear  part  of  the  equation, 

one  has  a  eigenvalue  problem  for  w,  in  which  the  following  determinant 
must  vanish  in  order  for  the  solution  to  exist. 


iw  +  vk2  ifk  /k2 

z 

ifk  iw  +  vk2 

z 

iRk  0 

x 


igkx/k2 

0 

iw  +  Dk2 


0 


(3.5) 


that  is,  (iw  +  vk2 )2 (iw  +  Dk2)  +  (fk  /k)2(iw  +  Dk2)  +  (Nk  /k)2(iw  +  vk2)-0. 

z  x 

This  equation  is  third  order  in  w  which  generally  yields  rather 
complicated  expressions  for  the  frequency  except  for  the  case,  D  «  v. 

We  will  briefly  indicate  below  the  effect  on  frequency  when  D  +  v  . 

For  the  remainder  of  the  derivation,  only  the  simple  case  D  ■  v  will  be 
considered.  However  the  procedure  used  below  for  obtaining  the  basis 
functions  remains  general. 


In  the  case  of  kinematic  viscosity  equal  to  diffusion,  the  solutions 
for  the  frequency  are 

-  ivk2 

u).  ■  w  +  ivk2  (3.6) 

2  o 

a),  — u)  +  ivk2 

3  o 

with  oj2  ■  (f2k2  +  N2k2)/k2  Thus  the  first  solution  is  the 

0  z  x  * 

diffusively  damped  zero  frequency  current  mode.  The  other  two  solutions 
are  frequencies  for  inertial  gravity  waves  with  dissipation.  When  D  ^  v, 
a  frequency  shift  from  the  natural  frequency  wq  and  additional  damping 
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are  to  be  expected.  We  briefly  Indicate  thle  by  considering 
asymptotic  behavior  of  u  near  e  «•  0,  with  e  -  (D  -  v)/v  .  Let 
Q  *  id)  +  vk2  .  Eq.  (3.5)  then  becomes  QJ  +  evk202  +  u>2Q  +  ev(fk  )2  -  0. 

O  2 

For  £«1,  this  can  be  solved  readily  by  perturbative  technique.  Letting 
Q  ■  Qq  +  2  +  ....  the  zeroth  order  solution  has  the  three  roots 

%1  *  0  ’  %2  "  i<>3o  ’  Q03  "  ‘  iU)o  ' 

To  the  first  order,  the  three  roots  are 

Q11  *  "  v(fkz)2/wo 

Q12  -  -  v(kxN)2/2w2  and  Q13  «  Q12 
To  the  second  order,  the  three  roots  are 
q2i-° 

Q22  "  i(vk2/2)(N/iJo)2(vk2/wo)[l  +  (3/4)(N/u)o)2(kx/k)2] 
and  Q23  -  -  Q22  • 

Thus  for  D  ^  v  and  e  «  1,  the  frequencies  are 

u.  *  iv[k2  +  tk2(f/u)  )2]  +  ... 

1  z  o 

w,  -  w  +  e2  (vk*/2)  (N/u>  )2  (vkz/u>  ) [1  +  (3/4)  (K/u  )*  (k  /k)2  ] 

“  o  a  o  o  ox 

+  iv[k2+  ek2 (N/u)  )z/2  +  ... 

X  o 

and  *  -  Real(ui2)  +  Imag(u>2) 

These  show  that  the  natural  frequency  tends  to  Increase  when  D  and  v  differ, 
and  the  dissipation  increases  or  decreases  slightly  depending  on  the 
sign  of  e. 

Now  returning  to  the  case  D  «  v  being  considered  for  the  derivation 
of  basis  functions,  the  dispersion  relations  are  (3.6).  The  eigenvector 
for  each  frequency  can  be  solved  for  and  these  are 
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xi.k '  (0-  11  h 


0.7*) 


*2,k  *  <1.  -'W  'RV"o) 


(3.7b> 


X3,k  "  <*•  W  «V“„>  * 


<3. 7c) 


where  the  elements  correspond  to  the  elements  in  the  vector  *  {iji.v.p'} 

A  A  A^ 

and  h,  a  ,  a  are  arbitrary  constants.  The  above  three  vectors  are 
linearly  independent.  Hence,  any  solution  of  ^jt»vit»Pjc^  is  gives  as  a 
linear  combination  of  the  three  vectors 


h  *  xi,k  e 


iu)lt  +  1  ,  el£U2t  +  I  .  e1U3t 


(3.8) 


At  this  point  we  may  generalize  to  arbitrary  t-dependence  by  denoting 
a^  *  a  exp(iu).t),  a^  *  a  er.pCiu^t),  h^  *  h  exp(iuijt)  , 

where  the  first  generalize  function  a^  may  be  thought  as  a  ’upward’ 
propagating  wave  mode  since  for  a  given  k  the  frequency  has  opposite 
sign  to  k.  Similarly.,  a ~  may  be  thought  of  a  ’downward’  propagating 
wave  mode,  and  h^  a  nonpropagating  current  mode.  Thus,  formally  we 
have  three  basis  functions  representing  three  elementary  me  les  of  motion 
in  a  stratified  rotating  fluid.  The  transformation  from  the  basis  to  the 
Fourier  coefficients  of  the  physical  variables  is  readily  obtained  from 
(3.8  )  .  in  matrix  form  the  transformation  is 


fVwo  -fkz/u,0  "gkx/fkz 


(3.9) 


[\\  [Rkx/wo  -Rkx/wo  1  J  l\\ 

It  follows  that  the  inverse  transformation  from  physical  variables  to 
the  basis  function  is 
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ak 
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(fk  )/(2k2w  ) 
z  o 

(&x)/(2k\) 

V 
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ak 

= 

]/2 

-(fk)/(2kV) 

z  u 

-(gkx)/(2k2£Do) 

\ 

\ 

0 

b 

-(fk2Rkx)/(lO30)S 

(fkz)2/(kwQ)2 

pk 

(3.10) 


Later  calculation  will  be  considerablly  simplified  if  all  the  above 
variables  are  nondimensionalized .  This  can  be  achieved  by  letting 
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\  - 


kR  + 
N  ak 


-  kR  - 
\  *  N  ak 


vl^i] 

jr  \  V  [l  \]  \  *  [r] 


(3.11) 
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ft 


t 


ft 


=  (fkz)/(kaJo) 


\  "  (Nkx)/(k^0) 


The  foregoing  transformation  matrices  in  terms  of  dimensionless  variables 
are 


V2  V2  V2 

\n  -0J2  -n^/2 

0  ‘  \  "k 


(3.12a) 


1/2k  1/lk  0 

“k 

\  "nk  “k 

k 


(3.12b) 
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C.  Some  identities: 


For  reference  purpose,  summarized  here  are  the  identities  resulting 
from  the  Fourier  transform. 

The  reality  condition  requires 


?.  -  4>  . 

k  -k 


* 

V.  -  V  . 
k  -k 


*' 

P  k  "  p-k 


(3.13) 


where  *  denotes  the  complex  conjugate. 

The  transformation  together  with  the  reality  condition  require 


K 


-k 


*-  + 
\  “  A 


-k 


(3.14) 


\  =  "H-V 


The  mean  energy  in  terms  of  Fourier  coefficients  is 


e,  -  i:  f.v.  +  v.v,  +  pr  p: 

k  k  k  k  k  k  k  k 


4-  *4-  —  *—  * 

2<\  \  +  \  \  >  +  »k  \ 


Finally,  for  the  dimensionless  parameters  we  note 


-k 


(3.13) 


-k 
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(3.16) 
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D.  Spectral  Equations  for  the  Basis  Functions 

Equations  (3.3)  for  the  Fourier  coefficients  can  now  be 
transformed  into  equations  for  the  basis  function.  This  is  done  simply 
by  substituting  M  for  and  then  multiplying  both  sides  of  (  3.3) 

by  M  ’  i’e‘’  M~1(W  -  A)  M  Y.  -  M_1  J.  .  N.  :  X  X 

Using  the  dimensionless  form  of  M  and  M  and  rewriting  the  W  matrix  as 


W  = 


“k  W‘t  1“kV*k 


‘VA 

‘■A\ 


0 

,2 


(3.17) 


by  applying  the  scaling  in  (3.11  ),  the  above  matrix  multiplication 
yields 

r. 


\  +  (iwk  + 

o 
o 


0  0 
2  > 


3C  +  (“ito^  +  y£)  0 

0 


< 

\ 

\ 

1 

-1 
(3.18) 


M-i  IN.  :X  X 

i  laa  ta 


where  we  have  scaled  w.  *  u)  /N  and  p,  «  vk2/NR2  .  The  nonlinear  terms 

k  o  k 

on  the  right  hand  side  of  (3.18)  are 

for  A*  :  1/2  I  .  i,  l  9.  ^  1(2,  H  ) 2..+  !1V  +  n.  p']  ,Q  ^ 

k.  £+a=k  k  £  k*£  £  £  £  a.  \  a  k  a  (3.19a) 

aT  :  1/2  E  .  2.2,  0.  ¥  [(£  /2  )2¥  2  -  ft.  -  T).p']  (3.19b) 

k.  £+£=k  k  £  k,£  £  £  5.  a.  k  k  £  k  £ 


Hk  : 


, e  .  2.2  9.  a*  i-a v  +  sip'] 

£+^=k  k  £  k,p  £  k  c£  li  £ 


(3.19c) 


where  9,  is  the  sine  of  the  angle  between  the  wave  vectors  k  and  £. 
_*£ 


Now  substituting  A+,  A  ,  and  H  for  H',V,p,using  the  transform 
relation  (3.12b)  and  expanding  the  nonlinear  terms,  we  have  for  the 
A^  equation  the  following, 
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3 A  *  «*v +  <i>  <  - 


(3.20) 


1/2  A-±  ViVt  (*1 +  V  1(VV  +  W  *  \VAa 

♦1‘W  ■  Va  •  Va1Aq 

-*W  Vi1  "a 

The  nonlinear  terms  can  be  simplified  somewhat  by  defining  a 
frequency  vector, 

^  -i\  ,  \}  (3.21) 

This  allows  rewriting 

Va  +  \na  *  V“a  ’  fk.i  :  \\  '  Vi  ‘  \,a 

which  are,  respectively,  the  cosine  and  sine  of  the  angle  between  the 

two  frequency  vectors.  Further  to  reduce  tie  number  of  terms,  we  introduce 

r  »  A+  +  A_  (3.22) 

£  £  £ 

which  is  the  sum  of  an  upward  going  and  a  downward  going  wave  of  wave- 
number  £  .  The  final  form  of  the  A^  equation  is 


3t*k  +  (iWk  +  *k 


(3.23) 


1/2  J.  .  a,  T  A  -  a  .  T  A  -  Y,.  „JH 

£+■£**  is, 23.  £  £  -ii. £9  £  £  L»££  £  £ 


where  the  interaction  coefficients  are 


a.  -  e.  (2  +  2.$.  ) 

£.££  iS.E  £  isii»£ 

Y,  »  2,  0,  9. 

iS*££  iS  iS»£  iS*£ 


(3.24a) 


(3.24b) 


Thus,  we  see  that  the  interactions  are  determined  by  the  ratio  of  the 
magnitudes  of  wave  vectors,  the  angle  between  the  wave  vectors  and  the 
angle  between  the  frequency  vectors;  Furthermore,  the  generation  of 
A+  waves  arises  from  both  wave-wave  interaction  and  wave-eddy  Interaction. 

In  a  similar  manner  the  A  wave  equation  can  be  obtained  by 
expanding  the  nonlinear  terms  and  the  result  is 


(3.27) 
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D.  Second  and  third  moment  spectral  equations 

In  the  foregoing  the  spectral  interaction  equations  for  A+,  A~,  and 
H  are  derived.  As  a  final  item,  we  will  write  down  the  governing 
equations  for  the  second  moment  of  these  three  functions  in  anticipation 
that  these  equations  will  be  useful  for  a  statistical  treatment  of 
rotating  stratified  wave-turbulence  interaction.  The  closeure  for  these 
equations  frequently  requires  the  consideration  of  the  third  moment 
equations.  For  completeness,  we  will  also  record  the  details  of  the 
third  moment  equations  at  the  end. 

The  three  basis  functions  A+,  A  ,  and  H  imply  six  different 
combinations  of  second  moments  among  them.  However,  we  will  record  only 
three  moments  here,  namely,  <A^  A*>,  <A^  Hm>,  <H^  Hm>,  since  the  moment 
equations  associated  with  the  A “  function  can  be  obtained  by  simply 
changing  the  sign  of  the  interaction  coefficients  as  will  be  indicated 
below.  The  angle  brackets  here  denote  ensemble  average.  Note  that  in 
the  case  of  homogeneity  which  is  often  assumed  in  statistical  study, 

•f 

the  number  of  moments  is  reduced  one  more  to  only  <Ak  V  and  <Hk  v- 

The  second  moment  equation  for  A.  A  is  obtained  by  forming  the 
■f  +  +  +  —  — 

products  A,  8  A  and  A  8  A,  .  The  two  products  are  then  added  and 
K  t  ID  m  t  K 

averaged.  The  same  procedure  applies  to  the  other  two  moment  equations. 
The  results  are  given  below 
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(3.28) 
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Now  note  that  the  spectral  equations  (3.23)  for  A  and  (3.25)  for  A 
differ  only  in  the  sign  of  their  coefficients  and  subscripts.  Thus  the 
above  second  moment  equation  for  A+  can  be  converted  to  one  for  A  by 
applying  this  simple  rule:  Change  the  superscript  of  the  A  function (s) 
in  the  time  dependent  term  to  negative  and  then  note  the  subscript 
associated  with  that  function.  In  the  remaining  terms,  change  the  sign 
of  the  frequency  associated  with  that  subscript,  the  sign  of  the 
interaction  coefficients  containing  that  subscript,  the  sign  of  the 
subscript  itself  in  the  coefficients,  and  the  sign  of  the  superscript  of 
that  particular  A  function  throughout.  This  completes  the  general 
second  moment  equations.  Among  the  second  moment  equations, of  particular 
interest  are  moment  equations  describing  the  spectral  energy  evolution. 
Eq. (3.15)  gives  the  expression  of  spectral  energy.  Since  A^  A^  •  A^  A_k 
the  above  equation  (3.28)  can  be  immediately  converted  to  one  describing 
the  spectral  energy  of  A^  mode  by  applying  the  just  stated  conversion 
rule.  The  left-hand  side  of  (3.28)  is  irreducible,  but  the  right-hand 
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side  is  simply, for  a£  A*  , 


for  \  K 


and  for 


w 


3t<Ak  ^k*  +2yk<Ak  A-k> 


3t<Ak  AV +  2i,k<Ak  Atk> 


-3t<Hk  V  -  2uk<ak  "-k" 


The  right-hand  side  of  (3.28)  to  (3.30)  contained  third  moment 
quantities.  We  conclude  this  report  by  listing  four  third  moment 
equations  for  different  combinations  of  A*  and  H  basis  functions. 

Moment  equations  involving  A  functions  again  can  be  deduced  by  changing 
signs  according  the  rule  given  earlier. 
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APPENDIX 


Flow  Research  Technical  Memo  No.  236 


To:  Jim  Riley,  Mike  Vleissman,  Ralph  Metcalfe  Date: 

From:  Greg  Holloway 

Subject:  "Basis  functions  for  analysis  of  stratified,  rotating 
turbulence" 


Date:  April  1981 


The  object  of  this  note  is  to  set  out  the  equations  of  motion  for 
an  unbounded,  incompressible  fluid  which  is  uniformly  rotating  and 
uniformly  stratified  in  the  mean.  Especially  important  is  that  the 
derivacijn  remain  fully  valid  for  finite  amplitude  up  to  the  Boussinesq 
approximation.  V.’e  require  such  a  confident  basis  before  proceeding 
toward  a  statistical  dynamical  theory  for  the  evolution  of  internal 
wave/stratified  turbulence  fields  in  oceans,  atmospheres  and  labora¬ 
tories.  Calculations  which  follow  are  "straightforward  but  tedious"  as 
is  appropriate  when  laying  a  foundation.  Sufficient  detail  is  given  so 
that  no  significant  algebraic  effort  has  been  left  to  the  reader's 
imagination. 

1 .  Derivation  of  Interaction  Equation:; 

The  instantaneous  fluid  state  may  consist  of  finite  amplitude 
departures  from  the  mean  state  which  we  seek  to  describe  in  a  non- 
perturba: ive,  complete  way.  In  fact,  we  will  make  the  Boussinesq 
approximation  on  the  assumption  that  fluid  accelerations  are  small 
compared  with  gravity  and  that  vertical  excursions  are  small  compared 
with  the  density  scale  height.  Then  we  begin  by  writing  "usual"  field 
equations  for  velocity  u,  density  deviation  p  and  pressure  p  . 


3^  u  +  f_  x  u  +  u  •  V  u  +  V  p  f-  £p  -  vv  u  *  0 


(l.a) 


+  u  •£()  -  Ru^  -  y  V  p  -  0 


(l.b) 


(l.c) 
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Here  c  is  departure  from  a  reference  density  p  (xO  *  1  -  Rx-  with  R 

2  o  j 

a  positive  constant  and  gR  ■  N  ,  the  squared  buoyancy  frequency  where 
we  have  nondimensionalized  density  by  the  mean  density  at  the  mean  depth 
which  is  taken  as  x^  •  0  .  We  take  f_  parallel  to  £  upward  along 

x^  .  v  is  viscosity  while  y  is  a  density  diffusivity  resulting,  say, 
from  thermal  conduction  in  a  thermally  stratified  fluid.  We  will  consider 
this  to  be  che  only  source  of  stratification  and  so  will  not  consider 
double  diffusive  phenomena. 

Art  important  aside:  The  present  formulation  is  quite  different 
from  a  usual  way  of  setting  up  the  internal  wave  interaction  equations. 
That  method  considers  a  field  of  fluid  "particle"  displacements  (jt,0)£() 
or  to  some  "equilibrium"  position  £  (x,t)  .  The  two  concepts  can  be 

confused,  i.e.,  equivalenced.  Diffusion  is  not  considered.  Interaction 
equations  are  then  obtained  as  an  infinite  power  series  in  the  "small" 
displacements  |  _£  |  .  Basis  functions  obtained  from  this  series  do  not 
appear  to  be  capable  of  a  complete  description  of  the  possible  fluid 
states.  1  feel  this  approach  can  only  be  useful  for  very  small  amplitude, 
non-dissipative  waves. 

Returning  to  the  present  derivation,  we  can  evaluate  pressure  from 
the  divergence  of  (l.a): 


cx . 
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~j 3x  fu£ 


ox . 
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u£  uj 


£  ~  J  3xj  3xj 

is  the  unit  alternating  tensor  and  we  sum  on  repeated  indices. 


P  + 
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®  3xn 


p  -  0 


(2) 


'j£m 


where 

Next  let  all  deviations  from  the  mean  state  be  constrained  to  be  periodic 
over  some  large  length  scale,  which  may  go  to  infinity.  Nondimensionalizing 
that  length  scale  to  a  volume  of  2tt  ,  transform  from  the  continuous  spatial 
variable  x  to  an  infinite,  discrete  set  of  wavevectors  It  on  which  we 
denote  Fourier  coefficients  of  u  ,  p  and  p  by 


(  u 


.  v  /  ,  \  ik*x 
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with  components  of  u  .  denoted  u 

~  £  j.k 

From  (2) 
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,  £+£■* 
where  k~  -  k*k  .  Substitute  this  into  the  transform  of  (l.a)  to  obtain 
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Combine  some  terras: 
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Likewise 
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where  for  symmetry  we  have  introduced 


M. .  -  ~  /k.  F,  +  k  F..  \ 

j4.n1  2  I  £  jm  m  ji  J 


Now  (l.a)  can  be  rewritten 
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while  (l.b)  and  (l.c)  are 
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Equations  (7.a,b)  constitute  four  equations  in  four  unknowns  ,  P^ 

having  eliminated  pressure.  We  may  still  use  incompressibility  (7.c)  to 
eliminate  one  component  of  velocity.  To  do  so,  define  a  set  of  unit  vectors 


for  each  k 


e£  -  k  x  e ' 3  /  |  k  *  e 3 


1  *  <  Ji  *  2  )  /  |  1  *(  k  *  eJ  ) 


e  *  V  x 

^k  -  3 


With  respect  to  these  unit  vectors,  define  new  velocity  components 


vk  “  -  k  -  k 
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*  k  * 
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The  horizontal  component  of  k^  is  k^  ■  k  *e^ 
to  denote  ratios  of  k^,  and  k  as 


It  will  be  convenient 


“ •  V  k  • e ■  V  k  >  »-lA 


(10) 


whence  uv  eliminate  u  k  by  incompressibility 


u  k  -  -  C  (k)  w  k 


(ID 


The  next  object  is  to  rewrite  (7.a,b)  in  terms  of  w  k  ,  v  ^  ,  p 
For  the  moment  we  will  not  worry  about  the  nonlinear  terms.  Then  we 
may  consider  a  single  wave  vector  k  and  so,  for  the  moment,  omit  the 
k  subscripts  as  being  implied,  u  k  can  be  written 

12  3 

u  *  v_e  -  0w£  +  we_ 
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We  will  need  £  ,  £  ,  £  in  the  original  Cartesian  coordinate  system. 

Their  components  are 
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where  s  *  1,2  and  k^  *  k^  +  k^  . 

Take  first  the  j  «  3  component  of  (7. a). 


Expand  the  term 
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Denoting  the  nonlinear  term  only  symbolically  as  iMuu  ,  the  j 
of  (7. a)  is 


3  component 
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To  obtain  an  equation  for  v  ,  contract  (7. a)  with  e*  .  This  is 
k. 
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Expanding  and  collecting,  obtain 
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So. 


(c*t  +  vk  )v  +  fow 


iCsj3  kj  (M  uu)s 


Finally,  from  (7,b) 


(14) 


(<*t  +  Yk  )p  -  Rw 
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where,  again,  we  write  the  nonlinear  terms  symbolically  and  will  worry 
about  this  later. 

Gathering  (13),  (14)  and  (15),  the  system  is 
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Now  we  can  check  that  for  linearized  waves  with  v  *  Y  *  0  ,  solutions 

proportional  to  e^o£  may  be  found  when  the  determinant  of  the  array 
in  (26)  vanishes.  This  yields 
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(17) 


as  well  as  a  solution  for  u-  ■  0  .  We  have  worked  very  hard  to  find 

o 

linearized  internal  waves,  but  at  least  this  tends  to  check  our  algebra! 
In  fact  (16)  is  valuable  because,  with  a  little  more  work,  we  could 
evaluate  completely  and  exactly  the  nonlinear  terms  on  the  right  side. 
This  is  in  contrast  to  recent  wave-wave  interaction  theories  for  which 
nonlinear  terms  are  given  by  an  infinite  power  series  which  cannot  be 
evaluated  exactly  for  any  nonzero  amplitude  of  the  wave  fields. 

A  quick  overview  of  where  we  are  and  where  we  are  going:  A 
minimal  statement  of  the  dynamics  given  by  (1)  requires  three  equations 
in  three  unknowns,  as  seen  in  (16).  If  more  unknowns  are  present, 
there  is  some  diagnostic  relation  among  them.  If  fewer  unknowns  are 
present,  the  description  is  incomplete,  as  recent  wave  interaction 
theories  may  be. 


For  unknowns  we  have  considered  w  ,  v  and  p  at  each  k  .  Call 
these  a  "primitive"  basis  ^  *  (w,  v,  p>k  .  In  terms  of  £  ,  (16) 
has  the  form 

(  i  3t  +  h  >  1= 

where  I  is  an  identity  and  L  and  K  are  linear  operators.  Other 
bases  are  possible,  in  particular  a  "wave"  basis  z^  »  (a+,  a”,  h>k  , 
which  will  satisfy  an  equation  like  that  for  jr  but  with  a  diagonal 
operator  in  the  nondissipative  case. 

Were  we  to  solve  for  the  time  evolution  of  all  ^  (or  z^)  ,  this 
would  give  us  the  history  in  detail  of  the  flow  field.  We  suppose  that 
is  neitr.er  practical  nor  especially  useful.  Rather  we  imagine  an  ensemble 
of  statistically  similar  flows  and  we  seek  to  describe  the  ensemble 
average  (denoted  <  •  >  )  moments  of  the  flow.  Suppose  there  is  no 
ensemble  mean  motion,  hence  <^>  “  0  •  Interest  then  turns  to  the 
second  moment  tensors 

h  •  %  V  °r  h  *  %  V 

Scn<_  of  the  components  of  include  density  variance,  <  pp*  >k  ; 

vertical  mass  flux.  Real  [  <  pw*  >k  J;  and  enstrophy, 

<(-  X  — )3>  k  "  K!  <VV*\  ’ 
while  total  kinetic  energy  is 

2 

KE  =  (1  +  o  )  <wv*>  +  <w*> 

_k  _k 

Gravi tat ional  potential  energy  is  related  to  <  pp*  >k  as 

PE  -  <  pp*  >k  N2  /  R2 

Reality  of  p  ,  w  and  v  in  physical  space  is  expressed  by 
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Consequently ,  e.g.  <pp*>^  *  <pp*>_^  . 

Next  let  us  explicitly  obtain  the  wave  bases  z.  ,  omitting  dissipation 

—  2 

for  the  moment,  i.e.,  v  *  y  ■  0  .  From  (16)  see  (3  +  wq)  ■  0.  So  obtain 

a+  and  a  as 


(9  +  iu  )  w 

t  o' 


2  — 

*  fu£v  -  gs  p  +  iuQ  w 


which  then  satisfy  (9  l  i-'  )  a”  ■  0  .  The  third  basis  component,  denoted 

to 

h^  ,  is  ncr.-propagating  and  can  be  seen  from  (16) 

3t  (v  +  IT  p)  "  0 


so  h  ■  v  -  ~  p  .  Together  the  transform  is 
K 


r  ,  r  2 " 

- 

‘  +  ■ 

-i  w  fu£  -gi 

w 

a 

o 

2 

iu^  fu£  -g-i 

V 

- 

a 

i 

Dl 

**-«  I  o<i 

H 

o 

_ 1 

p 

■> 

k 

h 

(18) 


such  that  the  linear,  inviscid  system  is  given  by  a  diagonal  operator 

0  0 


o  +  iuj 
t  o 


a  -  iuj  o 

t  o 


(19) 


Let  us  just  note  the  inverse  transform: 
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i/2 

-i/2 

0 

fc 

fa 

2  2 
N  a 

2_ 

o 

2u 

o 

w  2 
o 

R 

R 

Rf  up 

2,j 

o 

2u 

o 

2 

L) 

Ph> sically,  what  are  a  ,  a  and  h  ?  The  propagating  modes  a 
and  a  can  be  termed  the  "upward”  and  "downward”  propagating  inertial- 
gravity  waves.  In  terns  of  vertical  phase  propagation,  the  linearized, 
inviscid  waves  are  described  by 

ik,x,  i(±io  t  +  k-x*) 

±  ,  .  33  o  3  3 

a  (t)  e  «  e 

so  that  with  upward  we  have  a+  the  upward  going  wave  in  the  sense 

of  phase  speed.  Recall  for  these  waves  that  the  vertical  component  of 
group  steed  is  of  sign  opposite  the  phase  speed  so  that  energy,  say, 
associated  with  a+  is  downward  going. 

Perhaps  more  interesting  are  the  h  modes.  To  characterize  these 
consider  the  pure  case  when  a+  *  a  *  0  .  Then  from  (20),  or  from  the 
simultaneous  conditions  a  =  a  *  0  ,  we  see  that  pure  h  inodes  are 
given  by 

fov  =  gb  ,  w  ■  0  (21) 


i.e.,  the  horizontal  velocity  field  is  geostrophically  balanced  by 
variations  in  the  density  field.  When  we  expand  the  right  hand  side 
of  (16)  in  wave  bases  z  z  ,  then  if  we  retain  only  h  h  components, 

“2  ”SL  “E  ~a 

we  will  express  "quasi-geostrophic  turbulence".  However,  we  will  also 
see  the  transfer  from  quasi-geostrophic  motion  into  the  inertial-gravity 
waves.  When  inertial-gravi ty  waves  are  excited,  we  may  also  see  energy 
transfer  from  the  waves  to  the  quasi-geostrophic  field.  This  is  most 
tantalizing,  suggesting  a  possible  synthesis  among  quasi-geostrophic 
turbulence,  internal  waves  and  stratified,  three-dimensional  turbulence. 

While  it  is  generally  acceptable  to  refer  to  h-roodos  as  gcostrophic 
modes,  two  other  specific  interpretations  are  interesting.  Firstly,  when 

4 

f  *  0,  kj  4  0  then  h  =  v  while  v  drops  out  of  a*  .  Vertical 
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vorticity  (V  x  u)^  ^  ■  k^v  is  then  given  solely  by  h-modes  which  we 
would  now  be  inclined  to  cell  "vortex  nodes"  although  they  are  the  f-*0 
limit  of  geostrophic  moes.  Are  these  the  horizontal  pancake  eddies 
left  behind  in  wakes  after  more  three-dimensional  disturbances  subside  -- 
perhaps  by  radiation  in  the  a'  modes?  Do  nonlinearities  in  (16)  lead  to 
transfer  of  energy  from  h-modes  to  radiating  a-  modes  during  wake 
collapse? 

A  second  point  when  f  ■  0  is  the  case  k^  0  .  In  this  limit, 
i.e.,  vertical  k  ,  o  may  also  contribute  to  h  since  fo  is  undefined. 
Tnis  vertical  variation  of  horizontally  uniform  density  is  termed  "density 
f inestructure".  The  waves-f inestructure  interaction  is  fully  included  in 
the  present  analysis. 

Another  case  f  0  ,  k^  -*  0  ,  corresponding  to  density  f inestructure 
with  rotation  present  appears  to  be  troublesome  insofar  as  h  -*  ®  . 
Probably  it  would  have  been  more  appropriate  to  define  h  differently  as 
a  variable  h  for  present  purposes. 


h  »  ah  *  av  +  —  p  (22) 

recognising  that  the  dual  limit  f  ■»  o  ,  k^  o  is  special,  h  like 

a'  have  been  defined  so  that  z^  z  has  dimensions  of  energy. 

This  is  as  far  as  seems  useful  to  go  in  developing  and  discussing 

the  different  bases  for  the  problem.  What  is  important  is  that  we 

retain  a  complete  basis  for  describing  the  possible  fluid  motions. 

In  particular  we  do  not  discard  the  geostrophic  modes  (or  "vortex 

modes"  or  "density  f inestructure")  as  has  been  the  rule  with  previous 

internal  wave-wave  interaction  studies.  Rather  we  view  as  very  exciting 

♦ 

the  possible  exchanges  among  h  and  a~  modes  on  all  scales.  To  see 
these  we  turn  to  expanding  the  right  side  of  (16)  both  in  primitive 
bases  and  in  wave  bases  . 


—  1 1  A— 


Expansion  of  the  Nonlinear  Terms 


Firstly  in  primitive  bases  we  develop  the  right  side  of  (13): 


l3£m,k  ^  u£,£  Um»3. 


To  do  so,  first  note 


M  *  k 

333, k  *3 


Denoting  by  b  and  c  indices  i  3  , 


M3b3.k  '  M33b.k  '  2  k3  +  5  H  ^  ^  ) 


3  S,  ^ 


2  -  e2 

k  k 


K  (•  •  i  (  ¥=) 


(23. a) 


(23. b) 


-2k  S,  kc  k 


Using  these  we  express  some  terms 


,  u,  u,  *  k.u  w 
333,  k  3,2  3, a  3  £  a 


M3b3,k  %,£  U3,a 


\  \  (uk  'e  k>  (tbj3  pj  V£  ~  °£  Pb  V  W£  Ph 

K'v  6i>  (sj3  £  Vs  -  °E  k.  Vs> 


(23. c) 


(24. a) 


(24. b) 
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Likewise 


M33b,k  u3,£ 


K 


-  &  kn 


it 


v  w  -  o 

a  £  a 


%  \ 

—  w  w  I 

qh  aa/ 


Next 


V  it 

3bc,k  %,£ 


u 

c.a 


•  -6,  k-k  k  *(c,  .  ,  p.  v  -  p.  o  w  ) 
k  “T)  c  b j 3  £  M>  £  £ 

*(c  ,  q  v  -  qo  w  )p.  *  q~* 

cm3  £  c  £  fi  b 

ei 

'  Ya  k^  <W,'c  -  S,  Vb> 


+  W£  V£  kp^  bh  S>  pb  °p  ecm3  **■ 
pk 

+  W  kph  %  ^  kc  Ebj  3  PJ  qc  °£ 

6k 

-  w  v  t -  o  0  p,  q  k.  k 

£  fi  kPh  %  £  S.  b  c  *>  c 


Collecting  terms,  the  right  side  of  (13)  is 

L‘3£m,k  ^  Ui,£  Um,£ 

£4£-k 

•  -i  ^  IcjCkiB^Jv^  +  C2(k,£,£)v£ 

£+£-k  I 

+C3(k,£,£)V£  ♦  C3(k,£,£)Va  W£  j 


(24. C) 


(24. d) 


(25) 
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vile  re 
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(e^  *  k  x  £)  (e^  •  k  x  3) 
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°£  °a  S  kc  Pb  ’c 
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kph  % 
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(®  *  W  X  £) 


For  the  right  side  of  (14)  we  need 


k. 


-it 


s  j  3  k^  ^s£m,k  ^  U££  Um£ 

£+£“k 


Now  s  4  3  and  again  b  and  c  denote  indices  f  3.  Note  that 


M  -  k. 

s33,k  3 


(-7) 


‘  ksEk2 
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M  .  ,  .  s  M  . 
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soc.k  2  d  sc  c  sb  .2 


We  begir.  evaluating. 


sj3  Ms33,ku3pu3£ 
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Likewise 
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(27. a) 
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Finally  we  develop  as  above  and  find 
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Sj3 


M  .  ,  u.  u 

sbc.k  bp  cc^ 
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Collecting  terms,  the  right  side  of  (14)  is 
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*  -i  ]>  |cA(k,£,£)v£v^  +  C5(k.£.a)w£« 

£4£«k  I 

+  c6(k.£,a)v£wa  +  c6(k,a,£)vaw£  I 


where 
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The  right  side  of  (15)  is  more  easily  seen: 
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Tog,*ther,  (25),  (28)  and  (29)  give  the  complete  right  side  of  (16).  In 
matrix  notation,  that  right  side  may  be  written 

1  X  ‘  *2 

where  N  is  a  3  x  3  x  3  matrix  whose  elements  are  functions  of  k  ,  £ 
and  £  .  For  each  k  ,  £  and  £  there  are  27  elements  of  N  .  However, 
only  ten  are  different  from  zero.  Explicitly  these  are 


Nlll  =  C2(-’^*^) 
N112  * 

N12i  =  c3<k*^»£) 
N122  "  Cl(i^-2-^) 
N:n  * 

^2  12  * 

N221  “  c6(-*a-£) 


*332 
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In  like  fashion  was  seek  Co  develop  the  interaction  equations  among 
a+,  a~  and  h  .  The  straightforward  but  frightfully  tedious  method  is 
to  make  the  substitution  (20)  throughout  (16). 
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